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THE EFFECTS OF AN AUTOPILOT ON AIRPLANE RESPONSES 
TO TURBULENCE WITH EMPHASIS ON TAIL LOADS 

By Boyd Perry III 
Langley Research Center 

SUMMARY 

An analytical study has been made to assess the loads developed on the horizontal 
tail of an autopilot -controlled rigid airplane flying in one -dimensional atmospheric tur- 
bulence. The root -mean- square values of rigid-airframe responses and tail-load re- 
sponses were calculated at five flight conditions, and the behavior of these responses was 
observed in two autopilot modes: pitch -attitude -hold mode and altitude- control mode. It 
was found that pitch attitude and altitude can be controlled by the simple autopilot with 
acceptable or no increases in tail loads. 

INTRODUCTION 

The functions of most autopilots in transport airplanes range from stabilizing motions 
that have poor or no inherent stability to maintaining the airplane on a prescribed flight 
path. In the process, the automatic -control system must compensate for changes in vehi- 
cle and flight conditions and for atmospheric disturbances. A number of studies have been 
conducted on the responses of autopilot-controlled airplanes to drafts and turbulence; 
however, these investigations were directed primarily toward responses in airplane mo- 
tions and load factor (refs. 1 and 2). Presumably, designers' studies consider the tail 
loads associated with automatically controlled flight in turbulence, but practically no in- 
formation appears in the literature on the subject. 

An exploratory study, therefore, was made to reveal possible problem areas with 
regard to loads on the horizontal-tail surface arising from autopilot-controlled flight in 
turbulence. For this purpose, the airplane, autopilot, and turbulence are approximated by 
the simplest mathematical models that are felt to retain the elements which significantly 
influence the tail loads. The airplane is representative of the class of small corporate 
jet transports and is idealized as a rigid body in the longitudinal mode. The autopilot is 
approximated by the use of only first-order servosystem dynamics (ref. 3) and features 
pitch -attitude -hold and altitude -control modes through deflection of the elevator. The 
atmospheric turbulence is assumed to be one-dimensional for this analysis, and the usual 
assumptions of the characteristics for short Gaussian samples are made. The responses 
to turbulence include tail-load responses (in the form of root-mean-square values of shear 
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and bending moment at the root of the horizontal tail and torque about the horizontal -tail 
elastic axis) and airframe responses (root-mean- square values of pitch attitude, altitude 
perturbation, and center-of -gravity (c.g.) normal acceleration). In this study the various 
airframe and tail-load responses are calculated as functions of appropriate control-system 
gains for two centers of gravity, two cruise altitudes, and three autopilot servosystem 
lag conditions. 


SYMBOLS 


c.g. 

“t 


'L,o 


C Lt 


m 


c i 


c t (y) 

C 1 


normal acceleration of airplane center of gravity, positive up 
normal acceleration of horizontal tail, positive up 
span of horizontal tail 

airplane lift coefficient at reference flight condition 
lift-curve slope of horizontal tail, 3 C-^ ^ 3a t 

lift coefficient of horizontal tail 

elevator lift effectiveness, 3 ^3 S 

pitching-moment coefficient 

longitudinal -force coefficient 

plunge-force coefficient 

mean aerodynamic chord of wing 

mean aerodynamic chord of horizontal tail 

effective moment-arm coefficient for force distribution -f^(x) in appendix B 
chord of horizontal tail as function of span 

effective moment-arm coefficient for additional lift distribution due to angle 
of attack 
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Cg effective moment-arm coefficient for additional lift distribution due to elevator 

deflection 

Cg effective moment-arm coefficient for basic lift distribution due to elevator 

deflection 

c^ effective moment-arm coefficient for mass distribution 

f frequency, hertz 

f damped natural frequency, hertz 

f matrix of vertical-gust forces 

g 

g acceleration due to gravity 

CO ) frequency-response function of response R 

h altitude perturbation 



L horizontal-tail additional lift due to angle of attack 

&rr 


L horizontal -tail additional lift due to elevator deflection 

a S 

L, horizontal -tail basic lift due to elevator deflection 

b 5 
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K w 


** (y) 


^ (x) 


* a (y) 
a s 


(x) 


V y) 


-e^x) 

^ t (x) 

-e t (y) 

M 

M b (y) 


m 


m t (y) 

q 

q 


R 


s 


s 


t 


total lift on horizontal tail, L„ + L + L, 

a a a S D S 

chordwise horizontal-tail additional -lift distribution due to angle of attack 

spanwise horizontal-tail additional -lift distribution due to angle of attack 

chordwise horizontal-tail additional -lift distribution due to elevator deflection 

spanwise horizontal-tail additional-lift distribution due to elevator deflection 

chordwise horizontal-tail basic -lift distribution due to elevator deflection 

spanwise horizontal-tail basic-lift distribution due to elevator deflection 

horizontal-tail length, distance from airplane center of gravity to horizontal- 
tail aerodynamic center 

general chordwise force distribution 

total chordwise horizontal -tail lift distribution 

total spanwise horizontal-tail lift distribution 

coefficient matrix of equations of motion 

horizontal-tail bending moment at spanwise location y 

mass of airplane 

spanwise horizontal -tail mass distribution 

matrix of generalized coordinates in equations of motion 

nondimensional pitch rate 

general response to atmospheric turbulence 

wing area 

horizontal-tail area 
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Laplace variable 


torque about horizontal-tail elastic axis at spanwise location y 
servosystem characteristic time 
forward-speed perturbation 
airplane forward speed 

nondimensional forward-speed perturbation, u/u Q 

horizontal -tail shear at spanwise location y 

perturbation velocity along z-axis, positive down 

vertical gust velocity, positive up 

chordwise coordinate 

spanwise coordinate 

spanwise location 

airplane angle-of -attack perturbation, w/u G 
gust angle of attack, vf 0 j u Q 

horizontal -tail angle -of -attack perturbation (see appendix B) 
wing angle -of -attack perturbation (a = a ) 

elevator-deflection-angle perturbation, positive trailing edge down 
downwash angle at tail due to wing 
damping ratio 


pitch-angle perturbation 



6 ref 
X EA 
x .25c 

M 


R 






<I> W (<^) 

w g 


desired pitch-angle perturbation 
sweep angle of horizontal-tail elastic axis 
sweep angle of wing quarter-chord 


nondimensional mass, m/p S _ 

/ 2 

atmospheric mass density 
root -mean- square value of response R 
root-mean-square value of vertical gust velocity 
transport time lag, '6j l / u 0 
power spectrum for response R 
power spectrum for vertical gust velocity 
unsteady lift function 
circular frequency, 2vl 
circular damped natural frequency, 27rf Q 

Nondimensional stability derivatives are indicated by subscript notation as follows 


3C x 

C 

3C x 

Bu 

a 

Ba 

9C z 

c 

3C Z 

B u 


Ba 

3C z 




3C 


/dc \ 
3 — 

l2 V 


3q 
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BC 


m 


®u 


Bu 


C 


m* 

q 



c 


L 

a t 




dC 


m 


Ba 



B C 


m 





3S 


Subscripts: 

a c.g./ g 

a t 

EA 

h 



L t 

M b 


T 

u 

V 


normal acceleration of airplane c.g. 
normal acceleration of horizontal tail 
horizontal -tail elastic axis 
altitude perturbation 

horizontal-tail additional lift due to angle of attack 
horizontal-tail additional lift due to elevator deflection 
horizontal-tail basic lift due to elevator deflection 
total lift on horizontal tail 
horizontal -tail bending moment 
horizontal -tail torque 

nondimensional forward speed perturbation 

horizontal-tail shear 

airplane angle -of -attack perturbation 
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a. horizontal-tail angle -of-attack perturbation 

5 elevator angle -of -attack perturbation 

6 pitch-angle perturbation 

METHOD OF ANALYSIS 
Statistical Representation of Responses 

The method for obtaining the airframe and tail-load responses to random atmos- 
pheric turbulence is explained herein. Random-process theory is a technique for statis- 
tically analyzing output or response data for linear systems. It relates the input (atmos- 
pheric disturbance) power spectral density function to the output (airframe and tail load) 
power spectral density functions through the various frequency-response functions 
(refs. 4 and 5). The power spectral density function (power spectrum) contains all the 
statistical information describing a Gaussian process, including the root- mean-square 
(rms) value. 

The expression for the rms value of general response, R, per unit gust velocity is 

ff R 


where a R is the rms value of response R (airframe response or tail-load response), 
o w is the rms value of the vertical component of gust velocity, and <i> R (o;) is the power 
spectral density function of response R. The upper limit of integration in equation (1) 
has a finite value in actual practice and is approximately equal to 200 radians per second. 
Also, for a linear system 

* R («> = 0 W U)|h r (c)| 2 (2) 

s 

where H R (^) is the frequency-response function of response R to flight in a gust field 
which is sinusoidal in the direction of flight, has variable wavelengths, and is invariant in 
the spanwise direction. (See appendixes A and B for development of H r (cj) from the 
equations of motion.) The Dryden representation of the one-dimensional atmospheric- 
turbulence power spectrum is given by 


$ R (w)do 


11/2 


w. 


(1) 
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The Dryden spectrum is chosen over the more accurate von Karman spectrum because of 
its simpler mathematical form; for the purposes of an exploratory study, the Dryden 
spectrum provides sufficient accuracy. For this study, the integral scale length L was 
chosen to be 762 m. 

In conjunction with the equations of motion and the tail-load equations, which are 
discussed in appendixes A and B, equations (1), (2), and (3) are employed to obtain the 
rigid airframe and tail-load responses (in the form of rms values per unit rms gust 
velocity). 


Mathematical Models 

A rigid airplane, that is, one with rigid-airframe degrees of freedom but no elastic 
degrees of freedom, is assumed for this analysis. Its motion is described by the longi- 
tudinal degrees of freedom: forward speed, plunge, and pitch, featuring phugoid and short- 
period modes. The automatic -control system (autopilot) is idealized by assuming first- 
order servosystem dynamics. The idealized longitudinal autopilot consists of a pitch- 
attitude -hold mode and an altitude -control mode. Only one autopilot mode is operated at 
a time and no automatic trim device is considered. The equations of motion of the airplane 
autopilot system are basically the same as the classical equations of dynamic stability 
except that effects of the autopilot are included and appropriate gust forces are added. No 
unsteady-flow effects are included other than for transport time lags to account for the 
phase shift in the gust wave between wing and tail. 

Figure 1 contains block diagrams of both autopilot modes. The attitude -hold auto- 
pilot is a pitch-displacement-type autopilot with a pitch-rate feedback loop for damping. 
Pitch-displacement feedback gain YL Q , and pitch-rate feedback gain , are the control- 
system gains. The altitude -control autopilot contains an altitude-perturbation feedback 
loop with control-system gain K h (altitude -perturbation feedback gain). The "airplane- 
dynamics’' boxes in the block diagrams contain the longitudinal equations of motion, and 
the "atmospheric -turbulence" boxes contain the gust forces. The "elevator-servo" boxes 
contain the transfer function describing the first-order servosystem dynamics. The 
mathematical description of these "boxes" is found in appendix A. 
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Tail Loads 


The tail loads under investigation are shear and bending moment at the root of the 
horizontal tail and torque about the horizontal-tail elastic axis. The tail loads all consist 
of an aerodynamic component arising from airplane motions, gust velocities, and control 
deflections, and an inertial component arising from accelerations at the tail. Expressions 
for other responses, such as lift on the horizontal tail and accelerations at the tail (which 
are needed in defining tail loads), are listed in appendix B with the derivations of the tail- 
load equations. 


AIRPLANE AND FLIGHT CONDITIONS 

Airplane characteristics and flight conditions for this study are shown in tables I 
and II. Table I lists characteristics of the example airplane. The variations in cruise 
altitudes, centers -of -gravity, and servosystem dynamic conditions for each of the five 
flight conditions examined are listed in table II. Table III lists the airplane stability 
derivatives for each flight condition. 

RESULTS AND DISCUSSION 

Stability boundaries in terms of control- system gains and servosystem character- 
istic times, for both autopilot modes, are contained in table IV. The stability boundaries 
are determined by the familiar root-locus method, in which the exponentials describing 
the transport time lags are approximated by the first two terms of the exponential-series 
expansion. In the forced-response analysis, in which the transport time-lag exponentials 
are retained, the characteristic rapid increase in turbulence response with increasing 
autopilot gain identifies the approach of dynamic instability. The gains corresponding to 
the response peaks in the forced-response analysis are very nearly equal to the gains 
defining the stability boundaries determined by the root-locus method. Turbulence re- 
sponses are obtained for only those combinations of gains and characteristic times which 
result in a stable airplane-autopilot configuration. The turbulence responses are found 
by employing equation (1) after the frequency responses for G , h, M^(y), etc. have been 
obtained by solving the equations of motion and the tail-load equations. 

Attitude-Hold Mode 

The purpose of the attitude-hold autopilot is to maintain pitch attitude through pitch- 
displacement and pitch-rate feedback. The effects of varying K e and K- Q on airframe 
and tail-load responses and additionally on the characteristic oscillatory motions are inves- 
tigated for various flight conditions. The responses of interest in this autopilot mode are 
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pitch angle, c.g. normal acceleration, horizontal-tail bending moment, and horizontal-tail 
torque. Because of the similarity between the shear and bending- moment response curves, 
only bending-moment and torque response curves are presented. However, any discussion 
pertaining to bending moment also pertains to shear. 

The characteristic motions for zero servosystem lag time are found to be stable 
for all values of displacement and rate gains. For small values of Kg and Kg the 
characteristic motions consist of the oscillatory phugoid and short -period modes. As Kg 
becomes large, the phugoid mode becomes overdamped (£ > 1) and transforms into two 
modes: one with a slowly decaying predominantly airspeed variation and the other with a 
rapidly decaying predominantly plunge motion. The short -period mode remains under- 
damped (£ < 1) but is transformed into a mode with a predominantly pitching motion 
as K q becomes large. The damped natural frequency of this pitching mode varies in 
proportion to \Jk^ for large Kg, as shown in figure 2. 

The characteristic motions for nonzero servosystem lag times are not stable for 
all values of displacement and rate gains. The pitching mode becomes dynamically unstable 
at progressively lower values of displacement gain as the lag time is increased. The 
values of K d at instability are increased, however, by increasing the values of rate gain, 
as indicated in table IV. 

Figures 3 to 7 contain plots of rms values of airframe responses (denoted by sub- 
scripts 6 and a„ /g) and tail-load responses (denoted by subscripts IVL and T) 
against displacement gain for the various flight conditions in the attitude-hold mode. The 
responses are plotted on the ordinate and the pitch-displacement feedback gain, Kg, is 
plotted on the abscissa, which has a modified logarithmic scale broken near the left end to 
accommodate Kg = 0. 

The five curves on each set of axes in figure 3 correspond to five values of rate gain, 
Kg . The servosystem lag time equals zero for all curves. Pitch response decreases 
with increasing displacement gain, for all values of rate gain, because it is the controlled 
quantity. Large values of Kg provide the smallest pitch response but the largest c.g. 
normal acceleration and tail-load responses. The c.g. normal-acceleration response in- 
creases and then approaches a constant value about 30 percent higher than that for the 
basic airplane as Kg is increased. This trend reflects the loss of the airplane's natural 
tendency to pitch down into an up-gust as the pitching motion is suppressed. The c.g. 
normal -acceleration response at large Kg is essentially that of an airplane restricted to 
plunging motion only. The tail-torque response exhibits a characteristic similar to the 
c.g. normal acceleration except that the increase is approximately doubled. Analysis of 
the torque-load components indicates that the aerodynamic torque load due to elevator 
deflection (for this particular combination of elevator-chord to horizontal-tail chord ratio 
and elastic-axis location) is small compared to the other load components. The torque 
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load at large values of K 0 is governed by the plunging motion as is the c.g. normal 
acceleration. The tail bending-moment response increases slowly for Kg, < 10 and then 
rapidly thereafter. Analysis of the bending- moment load components indicates that the 
predominant contributor to bending moment is the elevator deflection. The rapid increase 
in the bending moment for Kg > 10 at smaller values of Kg, is due to (1) a reduction 
in the pitching-mode damping ratio caused by the transport lag in downwash, and (2) a 
spatial resonance effect from the phase shift in the gust wave between wing and tail. For 
a sinusoidal gust field having a half wavelength equal to the distance between wing and 
tail, the pitching moments from the wing and tail are additive. This phenomenon results 
in a maximum pitching motion for airspeeds that cause this wavelength to coincide with 
the natural wavelength of the poorly damped oscillatory pitch mode. This condition occurs, 
for example, at a value of Kg of about 300 and Kg of zero. Although not shown, the 
airspeed perturbation is not significantly excited by turbulence for Kg = 0 and the re- 
sponse decreases with increasing Kg . 

In general, the effect of adding rate feedback to the system is to increase the damping 
of the pitching motion and thus to reduce the magnitudes of all responses. 

Effect of changing the airplane center-of-gravity location. - Figure 4 contains re- 
sponses for flight conditions I and II. The effects of Kg on airframe and tail-load 
responses for two c.g. locations at constant altitude are investigated. For the data pre- 
sented in figure 4, Kg is held constant at a value of 10 and servosystem lag time is 
zero. As the c.g. is moved rearward (from 0.27c to 0.31c), the short-period damping 
ratio increases and the damped natural frequency of the short period decreases. The 
former tends to reduce pitch response while the latter tends to increase pitch response. 

The effect of increased damping predominates, however, and there is a slight reduction 
in pitch response for all values of Kg. The slight increase in c.g. normal-acceleration 
response occurs because the airframe's ability to "pitch into" gusts is reduced. The 
increase in bending-moment response with rearward c.g. shift is a result of the net in- 
crease in the aerodynamic component over the inertia component. The torque responses 
show that the rearward c.g. shift increases the torque about the tail elastic axis for small 
values of Kg and decreases the torque for large values of K^. 

Effect of changing cruise altitude .- Figure 5 contains responses for flight condi- 
tions I and m. The effects of Kg on airframe and tail-load responses for two cruise 
altitudes at a specific center-of-gravity location are investigated. For the data presented 
in figure 5, Kg is again held constant at a value of 10 and servosystem lag time is again 
zero. There are significant changes in airframe and tail-load responses due to an increase 
in cruise altitude. Over the entire range of Kg, reduced air density at the higher altitude 
results in an increased pitch response and a decreased c.g. normal-acceleration response. 


12 



Both the aerodynamic and the inertia components of tail loads decrease with an increase 
in altitude resulting in lower values of bending moment and torque at all values of displace- 
ment gain, K^. Figure 5 also illustrates that increasing Kg has less effect on changing 
the magnitude of the tail loads at the higher altitude than it does at the lower altitude. 

Effect of changing the servo-system characteristic time.- Figures 6 and 7 contain 
plots of responses against displacement gain for flight conditions I, IV, and V and show 
the effects of changing servosystem characteristic (or lag) time at one airplane c.g. 
position and one cruise altitude. Lag times of 0, 0,037, and 0.094 second are examined. 

In figure 6, K\ = 0 and in figure 7, = 10. The shapes of the response curves for non- 

zero lag times (flight conditions IV and V) follow the curves for zero lag time (flight 
condition I) very closely up to K q of approximately 1. At values of Kg between 1 and 
10 (depending on the values of and t the pitch-mode dynamic instability is en- 

countered. The responses then exhibit the characteristic rapid increase with increasing 
displacement gain and depart from the response curves for zero lag time. The effect of 
increasing servosystem lag times is to restrict the range of Kg available before the 
pitch mode becomes unstable. The restricted range of displacement gains also limits the 
amount of pitch control available from the autopilot. This limitation is not entirely detri- 
mental since the tail loads are smallest in this K g range. The effect of adding rate 
feedback to the system is beneficial from two standpoints: it increases the K & range 
(delays the onset of instability) and it reduces the magnitudes of all responses. 

The selection of the proper combination of displacement and rate gains requires the 

consideration of numerous criteria such as the effects of structural flexibility and unsteady 

>* 

flow, effects on ride quality, and handling qualities, as well as pitch -attitude disturbance 
and loads on the wing and the tail. These criteria are all dependent to some extent on the 
frequencies of the various natural modes of the airplane autopilot system. For the cases 
of nonzero servosystem lag times, the natural frequencies of the rigid airframe modes 
and servosystem modes are well separated from the frequencies which excite structural 
modes and from frequencies in which unsteady-flow effects are significant. These air- 
frame and servosystem frequencies also provide satisfactory ride and handling qualities 
for passengers and pilots for all combinations of gains within the stability boundary. For 
the highly idealized cases of zero servosystem lag time, however, the natural frequency 
of the pitch mode becomes very high with increasing displacement gain. In this higher 
frequency range, the effects of structural flexibility and unsteady flow must be considered. 
Flexibility effects may be significant for values of of 20 and higher. From figure 2, 

for equal to 20, the pitch-mode damped natural frequency is approximately 5 hertz, 

a frequency near the frequencies of the lowest structural modes of a small corporate jet 
transport. Also, in the higher frequency range there is a phase lag and attenuation of gust 
forces due to unsteady-flow effects. This phase lag and attenuation would result in turbu- 
lence responses lower than those presented for values of Kg of about 10 and higher. 
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It is anticipated that ride quality during flight in turbulence would be adversely affected 
by the use of values of Kg in the range from 5 to 10 unless large amounts of rate damp- 
ing, K*, could be applied. For this range of values of K e , the pitch- mode frequency 
falls in the range of about 3 to 5 hertz, which is a frequency range found to be very annoy- 
ing for human passengers. Handling qualities (how the airplane "feels" to a pilot) are also 
affected by the choice of gains. Generally, pilots prefer an airplane with a pitching-mode 
natural frequency that is in the vicinity of 1/2 to 1 hertz, together with a nearly critically 
damped motion. As indicated in figure 2, the associated value of K e may be about 1.0 
or less. The corresponding value of damping ratio is only 0.22, however. Adding K.-j to 
the system provides additional damping preferred by pilots and lowers the magnitudes of 
all responses as well. 


Altitude -Control Mode 

The purpose of the altitude-control autopilot is to maintain a constant altitude through 
altitude -perturbation feedback. The effects of varying on airframe and tail-load re- 
sponses are investigated for various flight conditions. The responses of interest in this 
autopilot mode are altitude perturbation, c.g. normal acceleration, horizontal-tail bending 
moment, and horizontal-tail torque. 

The characteristic motions for zero servosystem lag time become unstable at a 
finite value of K^. For small values of K h , the phugoid mode is essentially the same as 
that for the basic airplane but, as increases, the phugoid mode transforms into a 
slowly decaying long-period oscillatory mode with a strong pitching component and lesser 
amounts of airspeed and angle-of-attack components. It is this long-period phugoid-type 
mode which becomes dynamically unstable with increasing altitude-perturbation feedback 
gain. For all values of K h corresponding to a stable system, the short-period mode 
remains essentially unchanged. Also, for nonzero values of K^, a slowly decaying mode 
of motion appears that is predominantly a variation in airspeed. 

The characteristic motions for nonzero servosystem lag times become dynamically 
unstable at progressively lower values of as the lag time is increased. 

Figures 8 to 11 contain plots of rms values of airframe responses (denoted by sub- 
scripts h and 9 ) and tail-load responses (denoted by subscripts Mb and T) against 
altitude perturbation feedback gain in the altitude -control mode. The responses are plotted 
on the ordinate and the altitude -perturbation feedback gain, K^, on the abscissa. 

The effects of on airframe and tail-load responses for one c.g. position, one 
cruise altitude, and zero-lag time (flight condition I), are shown in figure 8. The value of 
altitude -perturbation response for = 0 ("basic airplane" condition) is theoretically 
infinitely large in the absence of the consideration of air-density effects on propulsion. 
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This infinite excursion in altitude is never realized, of course, but it does indicate that 
without altitude control a small rigid airplane tends to behave like a particle in the very 
low-frequency samples of atmospheric turbulence. The minimum value of altitude pertur- 
bation occurs at K h = 5.5 x 10“ 5 , and then the response increases rapidly as the critical 
value for instability is approached. The c.g. normal acceleration and tail-load responses 
are invariant with K h up to the onset of instability, at which point they too increase 
rapidly. The reason for the invariance of the acceleration and tail-load responses with 
gain at the small values of gain is the following: these responses are produced essentially 
by the response of the short -period mode which is not changed with changes in altitude- 
perturbation feedback gain. Thus, there is a range of in which the altitude is con- 
trolled effectively, the normal acceleration and tail loads remain unchanged, and an ade- 
quate margin from the critical value for instability is maintained. 

Effect of changing the airplane center-of-gravity location.- Figure 9 contains re- 
sponses for flight conditions I and II. The effects of on airframe and tail-load 

responses for two c.g. locations at constant altitude and zero servosystem lag time are 
investigated. Since the response of the short-period mode determines the magnitude of 
the responses (except the altitude perturbation) for the small values of K^, the effect 
of a rearward c.g. shift on the short period is of interest. As mentioned in the attitude- 
hold analysis, moving the c.g. rearward has the net effect of reducing pitch response in 
the short-period mode which tends to increase the c.g. normal acceleration and tail-load 
responses. This trend also appears in the altitude -control mode for the entire stable 
range of K h . In addition, moving the c.g. rearward decreases, by almost a factor of two, 
the critical value of K h for instability. At the rear c.g. position, then, only half the 
range is available for altitude control. Even with a reduced K^ range, however, the rms 
value of attitude perturbation reaches a minimum value approximately equal to that attained 
at the forward c.g. position. The altitude -control autopilot is effective at both c.g. positions 
without increases in c.g. normal acceleration or tail-load responses. 

Effect of changing cruise altitude .- Figure 10 contains responses for flight condi- 
tions I and IE. The effects of Kh on airframe and tail-load responses for two cruise 
altitudes at a specific c.g. location are investigated. Again, servosystem lag time is zero. 
Increasing altitude causes a reduction in air density which significantly reduces c.g. nor- 
mal acceleration and tail-load responses for all values of Kh. The value of Kh which 
drives the long-period mode unstable increases by nearly a factor of four with increasing 
altitude and, as a result, more gain is required to control altitude at the higher altitude. 

Unlike the effects of changing the c.g., for which a single gain (K^ = 3.4 x 10 -5 ) 
provides satisfactory altitude control and an adequate margin from the unstable condition, 
the effects of changing altitude do not permit the choice of a single gain. A gain which is 
satisfactory at the lower altitude is not very effective at the higher altitude, and a gain 
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that is effective at the higher altitude will result in a dynamic instability as altitude is 
reduced. A need for a that is programed as a function of altitude is indicated. 

Effect of changing the servosystem characteristic time .- Figure 11 contains plots 
of responses against altitude -perturbation feedback gain for flight conditions I, IV, and V 
and shows the effects of changing servosystem characteristic (or lag) time at one airplane 
c.g. position and one cruise altitude. Lag times of 0, 0.037, and 0.094 second are exam- 
ined. The primary effect of increasing the lag time is to reduce the value of which 
drives the long-period mode unstable. Except near the onset of instability, the shapes of 
the response curves for nonzero lag times (flight conditions IV and V) are identical to the 
shapes of the response curves for zero lag time (flight condition I). Even with the reduced 
range, increasing the lag time does not affect the extent to which altitude is maintained. 
The minimum values of altitude perturbation response for nonzero lag times are nearly 
equal to the minimum value of altitude perturbation for zero lag time. 

The effects of flexibility and unsteady flow are not significant in the altitude control 
mode due to the very low natural frequencies associated with the controlled quantity. Ride 
and handling qualities appear to be adversely affected only when the dynamic instability of 
the long-period mode is approached too closely. 

CONCLUSIONS 

A study has been made to assess the loads incurred on the horizontal tail of an 
autopilot-controlled rigid airplane flying in one-dimensional atmospheric turbulence. 

The root-mean-square values of rigid-airframe responses and tail-load responses were 
calculated for three servosystem representations, two airplane center- of- gravity posi- 
tions, and two cruise altitudes in both pitch -attitude -hold and altitude -control modes. The 
behavior of the responses was observed for variations in the three autopilot gains. 

Attitude -Hold Mode 

The following conclusions can be drawn from the results of the attitude-hold analysis: 

1. Increasing pitch-displacement feedback gain, K^, reduced the rms pitch 
response but increased c.g. normal-acceleration and tail-load responses. 

2. In general, the addition of pitch-rate feedback gain, , added damping to the 
system and reduced the rigid-airframe responses and tail-load responses over the range 
of K e . 

3. The shapes of the response curves and their relative magnitudes do not change 
with changes in c.g. position or changes in cruise altitude. 
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4. First-order time-lag representation of servosystem dynamics resulted in an 
unstable pitch mode with increasing K 0 . The instability limited the amount of K# use- 
ful for control. 

5. A combination of K„ and K- is available for all servosystem characteristic 

a u 

times which reduces pitch response with very little degradation in c.g. normal acceleration 
and tail-load responses relative to those of the "basic airplane" condition. 

Altitude -Control Mode 

The following conclusions can be drawn from the results of the altitude- control 
analysis: 

1. For all low-altitude conditions investigated, with zero servosystem lag time, 
a value of altitude -perturbation feedback gain, K^, exists which reduces altitude- 
perturbation response to a minimum without increasing c.g. normal acceleration or tail- 
load responses. For the high-altitude condition, with zero lag time, the tail-torque re- 
sponse increases about 25 percent at the gain corresponding to minimum altitude pertur- 
bation response. 

2. Increasing eventually results in the long-period mode becoming unstable. 
Decreasing altitude, moving the c.g. rearward, and increasing servosystem characteristic 
time all result in reduced K h range because of the long-period mode instability. 

3. Gain must be changed with altitude to avoid the long-period mode insta- 
bility and to maintain effectiveness. 

4. Changes in c.g. position and cruise altitude have greater effects on the shapes of 
the response curves than changes in servosystem dynamics. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 24, 1973. 
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APPENDIX A 


EQUATIONS OF MOTION 


The equations used in the present analysis are essentially the classical equations of 
dynamic stability available from the literature (refs. 6 and 7) with certain refinements 
made to describe the specific system under consideration. The refinements to the classi- 
cal equations are: gust forces are added, aerodynamic forces are applied in a manner to 
account more accurately for the lag in wing downwash and vertical turbulence component 
acting at the tail, and effects of the autopilot are included. The axes system in which the 
equations of motion are written is the stability axes system with origin at the airplane 
mass center. 

Expressed in matrix form, the equations of motion are: 


[M] {q} 



(Al) 


where M is the coefficient matrix of q, the matrix of generalized coordinates, 4 > is 
the unsteady lift function, and fw g is the matrix of vertical-gust forces. The elements 
of equation (Al) are as follows: 
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APPENDIX A - Continued 


The generalized coordinates are 



The unsteady -lift function is assumed equal to unity over the entire frequency range. 
The matrix of vertical-gust forces is given by 



Except for the quantity in boxes in equation (A2), the terms in the first three rows 
and columns are the classical terms found in the literature. The fourth row and column 
of equation (A2) contains the autopilot terms and will be discussed subsequently. The 
terms in equation (A4) are the gust forces written in terms of conventional stability deriv- 
atives, and it is assumed that the derivatives C Z( j, C m q, C z - , and C m ^ represent 
effects of the tail only. The quantity in boxes in equations (A2) and (A4) introduces the 
transport time lags of wing downwash and turbulence velocity at the tail. 

Aerodynamics 

In the classical form of the equations of motion, the transport time lags of vertical 
gust velocity at the tail and downwash at the tail are approximated by retaining only the 
first two terms of the series expansions of the time-lag exponentials. In the present 
analysis the exponentials are retained, and equations (A2) and (A4) will reduce to the 
classical form if the exponential, e _TS , is approximated by 1 - rs. No other unsteady- 
flow effects are employed in the present analysis. The unsteady lift function, <p , is set 
equal to one over the entire frequency range. 

Autopilot 

The autopilot consists of motion sensors and hardware to drive the elevator. The 
elevator is constrained to deflect in response to pitch angle and pitch rate in the attitude- 
hold mode and in response to perturbations in altitude in the altitude- control mode. Only 
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APPENDIX A - Continued 


one autopilot mode is operated at a time and no automatic trim device is considered. The 
elevator deflection is defined as follows: 

Attitude -hold mode 


5 


1 

st ch + 1 



+ K^s) e 


(A5) 


Altitude-control mode 


s 


1 

stch + * 


V h 


(A6) 


where 

h ~ ~ ( a - &) 


(A7) 


and where represents the first-order servosystem dynamics, and t C h is the 

servosystem characteristic (or lag) time. Physically t c h represents the time required 
for the servosystem transient response to reach a value equal to ^1 - of its steady- 
state value and is typically well below the period of the short-period nfode. As the expres- 
sion for b appears in the equations of motion (both modes contributing), would equal 
zeio in the attitude-hold mode and and would equal zero in the altitude -control 
mode. It is assumed that the hinge moment due to gusts is negligible compared to the 
other terms in the elevator equation. The last element in equation (A4) is, therefore, equal 
to zero. 

Frequency-response functions describe the responses of the airplane to a unit sinus- 
oidal gust of varying frequencies. Rigid- airframe frequency-response functions are 
obtained by solving the equations of motion for the generalized coordinates, q, and then 
substituting iw for s. From equation (Al) 



Solving for the generalized coordinates 
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APPENDIX A - Concluded 


and since 0=1 for this analysis 



Substituting io > for s in equation (A8) yields the frequency-response functions, and 
dividing equation (A8) by w g gives the frequency responses in per unit gust-velocity 
form, Hjy w , where R denotes airframe responses u, a , etc., and tail-load re- 


sponses 


Mjj, and T 



The frequency-response functions of c.g. normal acceleration and altitude perturbation 
per unit gust velocity are 

Ha c. g ./g/V i “ ) = (A10) 


V« M = 4 U °[V 


The frequency responses of a and h are linear combinations of the frequency 

c*g. 

responses of a and 8 . 

Equations (A9) to (All), in conjunction with equations (1) to (3) in the main body of 
the paper, provide the rms values of the rigid-airframe responses. 
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APPENDIX B 


TAIL LOADS 

The tail loads under consideration in the present analysis are the shear and bending 
moment at the root of the horizontal tail and the torque about the horizontal- tail elastic 
axis. All tail loads are comprised of an aerodynamic component and an inertial component, 
and the tail loads are calculated using simple cantilever -beam theory. The equations 
presented allow for tail-loads calculations at any spanwise location y = y . The tail loads 
at the horizontal-tail root are obtained by letting y = 0. 

The lift on the tail consists of contributions from angle of attack and elevator- 
deflection angle and provides the aerodynamic component of tail loads. The tail angle of 
attack 




S 0 + a 


g 


u 



(Bl) 


contains motion terms and gust terms, both with appropriate transport time lags. The 
expression for elevator deflection angle has already been given in equations (A5) and (A6). 
The expression for total lift on the horizontal tail is 

L t = - pu o ^ f ' C L a 
2 \ a t 

1 2 

which reduces to _ pu Q StC L a t for the autopilot-off condition. The lift distribution is 

2 “t 

assumed semielliptical across the span, and the chordwise distributions will be described 
subsequently. 

The mass of the horizontal tail times the acceleration of the horizontal tail provides 
the inertial component of tail loads. The normal acceleration at the horizontal tail is 



a f = a„ „ 
t c.g. 


- l h s 1 2 e 


(B3) 


which contains contributions from c.g. normal acceleration and pitching normal acceler- 
ation. The mass distribution is assumed trapezoidal along the span and is assumed to be 
concentrated at the 55 -percent-chord position. 
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APPENDIX B - Continued 
Shear 

The shear on the horizontal tail is obtained by integrating the spanwise aerodynamic 
force distribution and the spanwise inertia force distribution over the span of the hori- 
zontal tail. Positive shear is defined to be in the positive z-direction. The shear at any 
spanwise location, y, is given by 


Vfr) = 




* t (y) d y 


^ t (y)dy 



m t (y) dy • a t 


(B4) 


where l t (y) is the spanwise lift distribution and m^Cy) is the spanwise mass distribu- 
tion, and L t is the total lift on the tail and a t is the normal acceleration at the tail. 

The quotient of integrals in equation (B4) represents the fraction of the tail lift outboard 
of spanwise station y. When y = 0 the value of the fraction is 1 and the aerodynamic 
component of shear is then simply half the lift on the tail. The inertial component of shear, 
when y = 0, is the product of half the mass of the tail and the acceleration of the tail. 


Bending Moment 

The bending moment on the horizontal tail is obtained by integrating the spanwise 
shear distribution over the span, or equivalently, by integrating the products of the aero- 
dynamic and inertial force distributions and their respective moment arms over the span. 
Bending moment is defined to be positive for a positive shear acting at a positive y arm. 
The expression for tail bending moment at any spanwise location, y, is 


M b (y) - 


V 2 

't t (y)(y-y)dy 
y=y L t 

b.l 2 

\ 

b t /2 2 

Jy=y 

^ t (y)dy 


J y=o 



(B5) 


The quotient of integrals in equation (B5) represents the distance from spanwise station y 
to the centroid of that portion of the spanwise lift distribution outboard of station y. This 
distance is the moment arm of the aerodynamic component of bending moment. 


Torque 

The calculation of torque about the horizontal -tail elastic axis requires two sets 
of integrations: integrations across the chord and then integration across the span. 
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APPENDIX B - Continued 


The chordwise integrations serve to establish the effective "moment arms" of the various 
chordwise force distributions from their "centroids" to the horizontal -tail elastic axis. 
The "moment arms" are expressed as a percentage of the horizontal-tail chord, and this 
percentage of chord is assumed to be constant along the span. The spanwise integration 
of the spanwise force distributions times their respective "moment arms" then provides 
the torque about the elastic axis at any spanwise location y. 

In performing the chordwise integrations, it is necessary to know how the presence 
of the elevator hinge affects the calculations. It is assumed that the elevator hinge is full- 
span, nonswept, continuous, and frictionless and can therefore transmit no moments or 
torques. The following sketches will serve to illustrate how the "moment arms," which 
result from the chordwise integrations, are obtained: 



In sketch 1, 't^(x) is a general chordwise force distribution over the horizontal tail 
which produces a torque, T, about the elastic axis (labeled EA). In sketch 2 the force 
distribution has been integrated and replaced by two concentrated forces Lj and L 2 
at distances and x 2 from the elastic axis. These forces, and L 2 , are equal 

to the forces forward and aft of the hinge line (labeled HL), respectively, and produce the 
same torque about the elastic axis as the original force distribution produced (assuming a 
fixed hinge in both instances). The torque about the elastic axis due to Lj is simply 
XjLj. The torque about the elastic axis due to L 2 (with fixed hinge) is x 2 L 2 which is 
equivalent to (x 2 + x 2 )L 2< When the hinge is assumed frictionless (no longer fixed), the 
torque x^ L> 2 cannot be transmitted across the hinge line; therefore, the only contribution 
to torque about the elastic axis due to L2 is x 2 L2. The torque about the elastic axis 
due to both Lj and L2 is 


T = x 1 L 1 + x 2 L 2 


x 1 L 1 + x 2 L 2 
CjL 


c t L 


= C i c t L 
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APPENDIX B - Continued 

where L is the sum of and L 2 and equals the chordwise integral of t ^x), c. 
is the effective moment -arm coefficient such that when multiplied by c^, the mean 
aerodynamic chord of the horizontal tail, the product is the effective "moment arm" of 
distribution l . (x). 

For purposes of calculating the torque about the horizontal-tail elastic axis, the 
chordwise lift distributions (due to a ^ and s ) are further broken down by using the 
notation and method described in the appendix of reference 8. The chordwise lift distri- 
bution due to angle of attack is described by ^ a (x), the "additional" lift distribution 

a 

due to a^, This distribution is equal to the distribution over the horizontal tail when the 

elevator is not deflected. The chordwise lift distribution due to elevator deflection angle 

is described by l (x) and l . (x), the "additional" and "basic" incremental lift 

a s D s 

distributions due to S . The sum of -i „ (x) and (x) represents the increase in 

a 5 D S 

tail-lift distribution when the elevator is deflected. These chordwise lift distributions are 
illustrated in figure 12, where i represents -t a or -t a and l b represents -i b . 

The quantities L (x), L (x), and L h (x) represent the chordwise integrations of 

H. 3. o D ^ 

a o b 

l , -t and ■£■ . , respectively, and the sum of L , L_ , and L, is the total 

V a s V a a a s b s 

lift on the horizontal tail, L t< References 8, 9, and 10 provide the method for obtaining 

L , L , and L, as fractions of L. for the particular elevator-chord ratio used 
a ’ a, ’ b „ t 

a o o 

in this analysis. 

It is now possible to write the torque equation. Torque about the horizontal -tail 
elastic axis is defined to be positive counterclockwise when viewed from the left side of 
the airplane. The torque at any spanwise location, y, is given by 


T(y) = 


b t /2 

c l c t(y) ' e a a (y) dy L a j 
y=y a a °- 

-f- 

bx / 2 
' 1/ 

c 2' 

y=y 

a/ 2 2 

N 2 

■t, (yi <Jy 


II 

o 

p' 

O 

II 

>> 

-> 

A / 2 


c 3 c t< y > A s (y)dy L b 
J y=y s 

-+• 

A / 2 

1 c 

A / 2 2 

J - 
J y=y 

■f b (y) dy 


Vo 5 



t a (y)dy 
a s 


(B6) 
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APPENDIX B - Continued 


The products c. • c t (y), where i = 1 . . . 4, are the effective "moment arms” of the 
chordwise distributions, and they are expressed as fractions of the local horizontal-tail 
chord. The terms l (y), -C a (y), and l b ^ (y) are assumed to be semielliptical 

across the span. 

Frequency-response functions of the tail loads and associated responses are 
obtained in the same manner as the frequency-response functions in appendix A. Fre- 
quency responses for tail angle of attack, tail lift, and tail normal acceleration are ex- 
pressed as linear combinations of the rigid- airframe frequency responses. The tail-load 
frequency responses are then expressed as linear combinations of tail lift and tail normal 
acceleration. 

From equations (Bl) to (B3) 


H 

a 



(ic.j) 



de 

da 


H 


a. I w. 


(iai) + ia> 


g 


u. 


1 

H */w >> 

1 S u 0 



H 


L, /w^ 
t / g 




H 


t/ w g 


(iai) + C. 



(B8) 


H 3 /w M = g H a 


’ a t/ w g |4 “ ,| Slla c.g./g/ w g (iw) + ' t h- 


(B9) 


The frequency-response functions for the tail loads become 


b t/ 2 

f ^ t (y) dy 

J — xr 


H V/w = 
/ s 


y=y 


A / 2 

* t (y) dy 

J y=0 


b t /2 

I \/ w g <U ' , 'l - m t (y)dyH a t /w g (i " ) 


(BIO) 
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APPENDIX B - Concluded 


A / 2 


-MyKy-y) dy 


y=y 




y=o 


-t t (y)dy 


/ w (i") 

2 L t/ W g 


hi 2 

m t (y)(y-y)dy H a ; w (iw) 

♦'ir— ^ 


y=y 


(Bll) 


b t /2 

f c 1 c t (y)^ a (y)dy 
J ' ^ 1 

J y=y 1 


Brp /u? iO^) = 


r/wj 


A / 2 

y-0 


| -e a (y) d y 


— H t / (i^) 

2 L a/ w S 


b t /2 

f _ c 2 c t ( y> t a (y)dy 

y=y 


b t /2 


y=° 


h t 


L ag/ W - (ic ^ 


(y) dy 


A / 2 

c 3 c t (yk b (y) dy 

A=y 1 

b t /2 

| ‘tW d » 

J v=0 8 


I U \l w e 


M 


A / 2 


f c 4 c t (y)m t (y) dy H , (i, 

J y=y 1 g 


(B12) 


Equations (B7) to (B12), in conjunction with equations (1) to (3) in the main body of 
the paper, provide the rms values of the tail loads and associated responses. 
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TABLE L- AIRPLANE CHARACTERISTICS 


Mass of airplane, m, kg 7860 

Mass moment of inertia about pitch axis, I , kg m 27 600 

Wing area, S, m^ 31.8 

Mean aerodynamic chord of wing, c, m 2.55 

Sweep angle of wing quarter -chord, X 25c> deg 28.6 

Nondimensional tail length, 'f./c 2.6 

Horizontal-tail area, S^, m2 7.15 

Mean aerodynamic chord of horizontal tail, c^., m 1.47 

Sweep angle of horizontal-tail elastic axis, deg 0 

Elevator chord ratio 0.35 

de j da 0.566 


TABLE II.- FLIGHT CONDITIONS 
[Mach number was equal to 0.75 for all flight conditions] 


Flight 

condition 

Cruise 
altitude, m 

c.g. 

position, 
percent c 

Servosystem 
characteristic time, sec 

I 

6100 

27 

0 

n 

6100 

31 

0 

m 

12 200 

27 

0 

IV 

6100 

27 

.037 

V 

6100 

27 

i 

.094 
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TABLE in.- STABILITY DERIVATIVES 


Stability 

derivatives 

Flight conditions 

I, iv, V 

II 

in 

Cx u 

-0.0642 

-0.0642 

-0.0993 

C x 

.039 

.039 

.086 

a 




C L,o 

.133 

.133 

.332 

c s 

-.171 

-.171 

-.427 

C z- 

-2.46 

-2.41 

-2.46 

a 




c z 

-5.62 

-5.62 

-5.93 

a 




c z 

-4.35 

-4.25 

-4.35 

q 




% 

-.472 

-.472 

-.472 

C m u 

0 

0 

0 

C 

rn • 

-6.47 

-6.18 

-6.47 

a 





-.841 

-.616 

-.841 

% 

-11.44 

-10.92 

-11.44 

c m 

m c 

0 

-1.24 

-1.21 

-1.24 
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Figure 10.- The effects of altitude-perturbation feedback gain, 1%, on 
airframe and tail-load responses for two cruise altitudes. Flight 
conditions I and III (table II). 
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(a) Chordwise "additional" lift 
distribution for a thin airfoil 
with a plain trailing -edge 
control surface. 



(b) Chordwise "basic" lift dis- 
tribution for a thin airfoil 
with a plain trailing-edge 
control surface. 



x/c 


(c) Chordwise lift distribution for a 
thin airfoil with a plain trailing- 
edge control surface. 

Figure 12.- Chordwise lift distributions. 
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